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Abstract 

Photoplethysmography (PPG) devices are widely used for monitoring cardiovascular function. However, 
these devices require skin contact, which restrict their use to at-rest short-term monitoring using single¬ 
point measurements. Photoplethysmographic imaging (PPGI) has been recently proposed as a non-contact 
monitoring alternative by measuring blood pulse signals across a spatial region of interest. Existing systems 
operate in reflectance mode, of which many are limited to short-distance monitoring and are prone to temporal 
changes in ambient illumination. This paper is the first study to investigate the feasibility of long-distance 
non-contact cardiovascular monitoring at the supermeter level using transmittance PPGI. For this purpose, a 
novel PPGI system was designed at the hardware and software level using ambient correction via temporally 
coded illumination (TCI) and signal processing for PPGI signal extraction. Experimental results show that 
the processing steps yield a substantially more pulsatile PPGI signal than the raw acquired signal, resulting 
in statistically significant increases in correlation to ground-truth PPG in both short- (p G [< 0.0001,0.040]) 
and long-distance (p G [< 0.0001,0.056]) monitoring. The results support the hypothesis that long-distance 
heart rate monitoring is feasible using transmittance PPGI, allowing for new possibilities of monitoring 
cardiovascular function in a non-contact manner. 


Introduction 

Photoplethysmography (PPG) is a non-invasive light-based method that has been used since the 1930s for 
monitoring cardiovascular activity mm®. These devices have been used to assess cardiovascular factors such 
as blood oxygen saturation, heart rate, autonomic function, and peripheral vascular disease |2|- Conventional 
contact PPG devices are fastened to the skin using a spring-loaded clamp onto peripheral sites including the 
finger, ear, and toe [2]. Standard devices are comprised of either a single or multiple light-emitting diodes (LED) 
and a photodetector, and monitor cardiovascular activity by evaluating the change in light intensity due to 
fluctuations in blood volume. However, contact PPG devices are limited to monitoring a single individual per 
device, resulting in a linear increase in cost per individual monitored, and limiting the number of concurrent 
measurements. Furthermore, these devices produce single-point blood pulse measurements, which does not 
provide an indication for blood perfusion patterns which may be important to tissue viability [4j. Finally, their 
contact nature makes them unsuitable for long-term continuous monitoring due to comfortability and movement 
artifacts [2]. 

To address these issues, non-contact photoplethysmographic imaging (PPGI) systems have been proposed. 
The goal of a PPGI system is to extract the blood pulse signal without tissue contact. Although existing 
designs differ, they are primarily comprised of the same primary components as contact PPG: a light source 
(LED) and a light detector (camera). Such systems rely either on active [5, S [T, IS ; 9j Oil] or ambient tissue 
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illumination mmmmm Both active and ambient PPGI systems are sensitive to temporal changes in 
uncontrolled ambient illumination. The severity of this effect on the system’s ability to extract the subtle blood 
pulse signal is apparent with studies employing data collection in dark room settings 0MHO]. Normalizing 
ambient illumination changes using software has been proposed m , however this technique relies on a spectral 
estimation of the ambient illumination, which may fail in difficult lighting conditions. 

A potential, yet largely uninvestigated, advantage of PPGI is the ability to perform long-distance monitoring. 
Long-distance monitoring was defined as monitoring at the supermeter level. Long-distance monitoring is difficult 
or infeasible with contact PPG devices, since they are either attached via a cable to a monitor, must store 
the data on the device, or must transmit the data wirelessly, resulting in the need for a specialized network 
infrastructure. Still, one device can only monitor one individual. In contrast, PPGI systems may acquire long¬ 
distance measurements of multiple individuals, provided that the individuals are in the camera’s field of view. 
However, long-distance monitoring becomes challenging with existing systems. Many existing PPGI systems 
operate in reflectance mode, where the camera and illumination are positioned on the same side of the tissue 
under investigation. In addition to the weakened reflectance signal by the highly scattering nature of skin mna, 
reflectance PPGI systems using divergent light sources are sensitive to a strong illumination intensity decrease. 
These factors render long-distance monitoring a challenging problem. These constraints may be alleviated using 
a transmittance PPGI system, where the LED is positioned close to a thin anatomical location (e.g., finger, toe). 
To the best of the authors’ knowledge, this paper is the first to investigate long-distance transmittance PPGI 
monitoring. 

This paper presents a pilot study to assess the feasibility of long-distance cardiovascular monitoring using 
transmittance PPGI. For this purpose, a novel non-contact transmittance PPGI system is proposed which is 
able to monitor cardiovascular activity remotely by correcting for ambient lighting fluctuations and extracting 
the subtle blood pulse signal using signal processing tools. This system comprises a 100 fps camera, a high- 
powered LED (655 nm), a microcontroller to synchronize frame captures and illumination, and image and signal 
processing software on a computer. Figure [l] shows a pictoral representation of this system, and Figure [2] shows 
the processing steps. Temporally coded illumination (TCI) is introduced to remove ambient lighting artifacts 
at the acquisition level, thus correcting the data prior to further processing, similar to that introduced in our 
previous work [IS]. The proposed PPGI system used for this study extends substantially beyond our previously 
reported work by incorporating a comprehensive set of signal processing steps, a larger testing sample size, and 
more rigorous statistical and qualitative analysis. The signal processing is particularly important for this study, 
as long-distance cardiovascular monitoring scenarios result in acquired signals with low signal-to-noise ratios 
(SNR). 


Results 

Experimental Setup 

For this pilot study, five of the authors’ data (age 28.2 ± 5.3) were collected using the proposed system. A 3:1 
TCI code (alternating 3 frames with LED on, followed by 1 frame with LED off) was used with 100 fps frame 
rate and 8 ms exposure time. The aperture was manually adjusted to limit pixel saturation. The EasyPulse 
PPG device [TB] was used to collect the ground-truth PPG signal for validation. 

Two experiments were performed to assess the feasibility of long-distance monitoring. In Experiment 1, the 
camera and LED were separated by 20 cm, serving as short-distance base case validation. The participants 
were asked to position their fingers between the LED and camera so that their fingers covered the beam of the 
LED. In Experiment 2, the camera and LED were separated by 1.5 m (“long-distance”), and the participants 
were asked to position their fingers at approximately 10 cm from the LED. Figure [3] shows example images of 
both experiments before and after ambient correction. For each experiment, a 10 s window was chosen that 
yielded a clean ground-truth PPG signal for validation. The normalized power spectral density (PSD) was 
computed for spectral analysis to demonstrate each signal’s dominant frequency components. In an ideal signal, 
the fundamental heart rate should be the dominant frequency. Furthermore, to assess temporal signal fidelity, 
the Pearson’s linear correlation coefficient p was computed between between PPGI and PPG signals. This metric 
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is offset- and scale-invariant, suitable for the problem of comparing unit-less PPG signals: 


p (x,Y) = ^-e[- 1,1] (i) 

a x cry 

where X,Y are the finger PPG and PPGI signal, ctxy is the covariance between the two signals, and crx, cry 
are the standard deviations of the same variables. This generates a value p £ [—1,1], where 1 is perfect linear 
correlation, —1 is a perfect negative linear correlation, and 0 is uncorrelated. A perfectly reconstructed PPGI 
signal would yield p = 1, indicating that the PPGI signal is a linear multiple of the finger PPG signal. 

The PPGI signals were statistically analysed using a novel physiologically-motivated f-test to quantify the 
degree to which the proposed system processing improved the signal. For each participant, each heart beat 
instance (i.e., diastole-to-diastole) was semi-automatically labeled using a custom gradient descent approach, p 
was computed for each such heart beat instance, yielding n correlation values. A two-tailed paired-sample f-test 
(a = 0.05) was conducted using the correlation values for each heart beat instance, yielding a p-value indicating 
the probability of a zero-mean difference between signal correlation. 

Experiment 1: Short-Distance 

The goal of Experiment 1 was to validate the base case short-distance setup required prior to long-distance 
measurement. Figure [4] shows the effect of processing the PPGI signal for three participants. The processed PPGI 
signal contained amplified pulsatility over the unprocessed (“raw”) PPGI. Pulsatile signals were extracted for all 
five participants (“PI” through “P5”), resulting in discernible systolic peaks of each blood pulse. The processed 
PPGI subdued much of the noise present in the raw PPGI signal. The high-frequency information not expected 
in the naturally smooth blood pulse waveforms was denoised, yielding smooth blood pulses. Furthermore, non¬ 
linear trends were corrected, yielding stable PPGI signals. These non-linear trends resulted in low correlation 
values between the raw PPGI and the ground-truth PPG (e.g., at 4 s for PI). The processed PPGI stabilized 
the signal, yielding larger correlation values at those time windows. The heart rate was discernible in both raw 
and processed PPGI as a peak frequency power coincident with the PPG peak. However, although all PPGI 
signals showed strong spectral similarity to the PPG signal in the PSD, the processed PPGI signal exhibited a 
stronger peak at the heart rate, thus subduing extraneous frequencies. 

Table [Tj summarizes the statistical results for this experiment. The correlation values across whole-signal 
comparisons were more than doubled for each participant from unprocessed (p = 0.24 ± 0.10) to processed 
{p = 0.72 ±0.19) due to recovered stability and pulsatility. This indicates a large correction across the entire sig¬ 
nal, demonstrating an overall improvement in signal fidelity. Statistical analysis using the heart beat instance cor¬ 
relation values showed a statistically significant increase in correlation for each participant (p £ [< 0.0001,0.040]). 
Thus, as a result of processing, the system was able to improve the PPGI signal such that there was a significant 
increase in pulsatility beat-to-beat, making inter-beat analysis more feasible. 

Experiment 2: Long-Distance 

The goal of Experiment 2 was to validate long-distance PPGI measurement for long-distance monitoring. Figure[5] 
shows the effect of processing the PPGI signal for three participants. The raw PPGI signals were noisy, and 
pulsatility was not easily discernible in the signals of PI and P3. This was emphasized by the lack of a clear 
peak in the PSD. Processing recovered the pulsatility for these participants, and yielded a less noisy and more 
stable signal. The processed PPGI signals yielded a higher degree of pulsatility, although the signals were noisier 
than the short-distance measurements (this is discussed later). This pulsatility was reflected in the PSDs, where 
there was a clear peak coincident with the heart rate for each participant. However, in cases such as P2 and P3, 
distinct blood pulses were not easily observed, making beat-to-beat timing analysis difficult at long-distances. 
The correlation plots showed a general increase across each participant with processing, demonstrating the 
increased pulsatility. However, there were some areas that were not improved (e.g., P3 at 2-3 s), largely due to 
unsuppressed noise likely due to amplified movement artifacts (discussed later). This high-frequency noise did 
not affect the PSD within the physiologically valid heart rate range. 

Table [2] summarizes the statistical results for this experiment. As in Experiment 1, the correlation values 
for whole-signal results more than doubled from unprocessed ( p = 0.14 ± 0.05) to processed (p = 0.55 ± 0.12), 
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indicating recovered pulsatility and stability for an overall increase in signal fidelity. Statistical analysis using the 
heart beat instance correlation values showed a statistically significant increase in correlation for each participant 
except P3 (p G [< 0.0001,0.056]). However, P3 exhibited substantial improvement upon processing (p = 0.056). 
This emphasizes the increased beat-to-beat recovery of pulsatile components. Overall, correlation values were 
lower than for short-distance monitoring, but significance (or near-significance) was achieved. 


Discussion 


Strong results were obtained for short-distance measurement results, producing visually pulsatile signals and 
accurately extracting heart rate. Statistical significance was achieved (p G [< 0.0001,0.040]), exemplifying the 
recovery of pulsatile PPGI signals from a transmittance signal acquired in a non-contact manner. The processed 
PPGI signal yielded smooth pulsatile components, making temporal beat-to-beat analysis feasible. For example, 
heart rate variability can be used as a measurement for autonomic nervous system activity |201I21] . The variability 
in correlation values across participants may be due to factors such as movement artifacts, and differences in 
scattering and absorptive media across tissues, such as melanin (skin tone) and adipose tissue |22] . Gender 
differences were not analysed, since the participants were all male. Overall, the results validate short-distance 
cardiovascular monitoring using PPGI. 

Long-distance monitoring yielded lower correlation values to ground-truth PPG than those for short-distance 
monitoring. Due to the diffuse nature of light traveling through highly scattering tissue [23] , the sensor irradiance 
is reduced, yielding weaker signals at greater distances. Furthermore, the spatial resolution of a pixel increases 
with distance, and thus small movement artifacts were amplified. This is apparent in the results for participant 
P3 (see Figure [5]). The participant’s fingers were shaky, yielding high-frequency noise which the system was 
unable to fully remove. During long-distance monitoring, statistical significance was achieved for heart beat 
instance correlation values in all but one participant (P3, p = 0.056), indicating the recovery of a pulsatile PPGI 
signal. 

Although pulsatility was extracted in some but not all cases, the fundamental heart rate was identifiable 
across all participants. This indicates that beat-to-beat timing analysis becomes difficult in the time domain, 
but heart rate is identifiable during long-distance PPGI monitoring. These results indicate the feasibility of long¬ 
distance heart rate monitoring, but further work must be done to demonstrate beat-to-beat temporal analysis 
in long-distance monitoring. Such an imaging system could benefit settings in which the individual exhibits low 
motion, such as cardiovascular monitoring during surgical procedures or sleep studies. In both settings, having 
as little contact on the individual’s body to limit the amount of obstruction is preferable. Such an imaging 
system could provide physiologically relevant information in a non-contact manner. Additionally, the imaging 
system could monitor multiple individuals within the camera’s field of view by analysing each individual’s signal 
independently. The cost per individual decreases with a greater number of individuals, since the setup requires 
a single camera, and one LED per individual. 


Methods 

To investigate the feasibility of long-distance PPGI monitoring, a non-contact PPGI system is proposed whose 
design integrates hardware (active illumination, a camera, and synchronization electronics) and software (signal 
processing, timing signals). Figure [l] depicts a graphical overview of the proposed PPGI system. A participant’s 
fingers were placed between a high-power LED (Philips LUXEON Rebel, 655 nm, 580 mW) and a camera (Point- 
Grey Flea3, 100 fps), enabling long-distance monitoring. The camera frames and LED illumination patterns are 
synchronized by a microcontroller (MCU) and custom electronics. Temporally coded illumination (TCI) is pro¬ 
posed for controlling this illumination pattern. Figure [2] shows the stages of processing to yield the final PPGI 
signal. The following subsections provide detailed description of the processing steps. 
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Light Transport Model 

A temporal extension of the Beer-Lambert law (BLL) was used. The BLL shows that light is exponentially 
attenuates as it passes through a homogeneous light-absorbing medium: 

T=^= e~ £lc (2) 

hi 

where T is transmittance, /o is incident light, I is transmitted light, e is the tabulated molar extinction coefficient 
for the light-absorbing medium (e.g., oxyhemoglobin), l is the photon path length, and c is the concentration 
of the medium. Noting that path length for oxyhemoglobin is the dominant temporally changing parameter in 
heart beat analysis [3], the temporal transmittance can be written as a function of time t: 

T(t) = ^P~ = e- £ - m - c (3) 

Iq 

The standard measurement signal for PPG analysis is absorbance, which is related to transmittance as A = 
— log(T). Then, expressing this as temporal absorbance measurement yields: 

A{t) = — log (I(t)) (4) 

This model assumes constant ambient conditions. However, I(t) may be affected by temporal changes in 
ambient illumination, potentially corrupting the subtle blood pulse signal. Thus, in the next section, an ambient 
correction method is proposed. 


Temporally Coded Illumination (TCI) 

A TCI process is proposed to remove the effect of ambient light changes in the data, thus ensuring that the 
reflectance signal is due solely to the controlled active illumination (e.g., LED). A dual-mode temporal code 
was defined by which frames contain irradiance due to either ambient illumination or a combination of ambient 
and active illumination. By coding the illumination synchronously with frame captures, ambient measurements 
can be performed according to the temporal code, and subtracted from the actively illuminated frames. For 
example, a temporal code of 3:1 defines a TCI sequence that measures one ambient frame after every three 
actively illuminated frames. Figure [I] shows a graphical example of 3:1 TCI. 

Mathematically, given an intensity image I(t) at time t , according to the TCI sequence, I(t) can either be: 


m = 


Iamb (^) > 

Iamb{t') T loot (^) 7 


LED off 
LED on 


(5) 


where I am b(t) and / act (f) are the irradiance due to ambient illumination and active illumination, respectively. 
Then, given a recently acquired ambient frame I am b(t— At), by modeling light reflectance as a linear combination 
of ambient and active lighting, the ambient correction can be applied to the current frame: 


I*(t) = I(t)-I amb (t-At) 


( 6 ) 


With a fast frame rate, At becomes negligible, yielding illumination due solely to controlled active illumination: 


lim I*(t) = I a ct(t) (7) 

Ai-vO 

Since the system camera operates at a frame rate of 100 fps, At approximately satisfies this condition (At = 
10 ms). Substituting ([ 6 ]) into @ yields the ambient-corrected absorbance measurement: 

A(t) = — log (/*(t)) ( 8 ) 

Once the ambient frame I am b(t ) is captured, the previously acquired active frame J(t— At) was used in its place 
to ensure constant frame rate for frequency analysis. 

TCI was implemented at the hardware level, and requires control of an active light source. Thus an active 
illumination system was chosen. Camera and illumination synchronization was performed using custom circuitry 
and microcontroller code. Figure [3] shows the results of ambient correction via TCI on a single frame. Media 1 
shows ambient correction across a video with varying lighting. 
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Signal Processing 

Capturing data using the given light transport model and TCI yields an ambient-corrected PPGI signal. However, 
this signal may still be affected by various sources of noise. This section proposes a signal processing pipeline 
for extracting a stable pulsatile PPGI signal. First, a denoising step is proposed, which reduces the noise based 
on camera and process characteristics. Second, a detrending step is proposed, which removes trends in intensity 
due to movement. 

Denoising 

The goal of signal denoising is to recover the smooth blood pulse waveform from a noisy signal. Various sources of 
noise may corrupt the signal, such as light shot noise, the inherently stochastic nature of light-tissue interaction, 
due to measurement noise due to camera read noise. Kalman filtering 1241 models observed measurements as a 
signal containing two general types of noise: process noise and measurement noise. This models the system well, 
considering the light-tissue interaction as process noise and sensor noise as measurement noise. It is a recursive 
algorithm that can operate in real-time, as denoising at time t depends only on the state and uncertainty matrix 
at time t — 1. Since blood pulse waveforms are inherently smooth in nature, the system state is defined using 
Newton’s laws of motion using the raw camera intensity values: 



r/fei 
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where, at time point k, 4 is pixel intensity, Ik is the rate of pixel variation, and Ik is acceleration of pixel 
intensity. The measurement noise was assessed by imaging a dark field and quantifying the standard deviation 
of pixel intensities a m across the sensor. During the experiments, a m = 1.02 for 8 bit pixel values and 8 ms 
integration time. Denoising was performed on the raw camera signal I(t ) rather than the computed absorbance 
signal A(t) since a m was learned using the untransformed raw camera intensity values. 

Process noise is more difficult to model due to the complex and stochastic interaction of light with tissue. 
Thus, a grid search optimization was performed to assess the optimal process noise variance value a p using 
8 ms exposure time. Using a fixed measurement noise model and a training data set, a logarithmic grid search 
was performed over various values for <j p . The parameter value that yielded the highest correlation to the 
ground-truth PPG signal was chosen, resulting in a p = 1.84 x 10~ 2 for 8 bit pixels. This denoising process was 
then performed on the acquired signal prior to any further processing. Then, the denoised intensity signal was 
extracted using an observation matrix: 

4 = Hx k ( 10 ) 

where H = [1 0] T . Plugging this into Q yields the denoised absorbance signal: 

A(t) = -log(I(t)) (11) 


Detrending 

Fluctuations in lighting and measurement or process noise should be corrected using the denoising algorithm. 
However, these corrections assume constant incident illumination, which is an invalid assumption during move¬ 
ment. The goal of this detrending step is to remove slow oscillations in the signal, yielding a PPGI signal with 
constant offset. It was assumed that movements are relatively smooth over time. A detrending algorithm was 
used on the denoised absorbance signal A(t), in which the model assumes a smoothness prior |25j . In particular, 
the observed signal A(t) is modeled as the linear combination of the “true” absorbance signal A true {t) and a 
temporal trend A tre nd{t). 

A(t) = Atr Ue (t) + A tre nd{t ) (12) 

Given that A(t) is measured, A true {t) can be solved by estimating A tren< i{t) assuming a linear model, subtracting 
it from A(t), and solving this using regularized least squares, which provides the following estimate of the true 
detrended signal: 

A true = (l~(I + X 2 D^D 2 )~ 1 )A (13) 
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where I is the identity matrix, A is a relative weighting term, and D 2 is the discrete approximation matrix of 
the second derivative. Figure [g] shows the recovery of a stable signal from a signal corrupted by movement using 
this detrending method. 
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Figure 1: The proposed non-contact PPGI system comprises an LED (655 nm), a CMOS camera (100 fps), and 
a microcontroller (MCU) for synchronization and TCI implementation. The MCU timer dictates the TCI code, 
triggering the camera to capture frames, which are transmitted and processed on the computer. 
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Figure 2: Overview of the processing steps for the proposed PPGI system. Upon capturing frames, ambient 
correction using TCI is performed, which removes temporal changes in ambient illumination from the frames. 
This is followed by denoising to remove camera sensor noise and process noise from the light-tissue interaction, 
and detrending to provide a stable PPGI signal. The resulting signal is an extracted stable PPGI signal. 



(a) Experiment 1, unprocessed 


(b) Experiment 1, ambient corrected 



(c) Experiment 2, unprocessed 


(d) Experiment 2, ambient corrected 


Figure 3: Example images for Experiment 1 (short-distance monitoring) and Experiment 2 (long-distance mon¬ 
itoring). The unprocessed frames (first column) contain uncontrolled ambient illumination (windows, overhead 
lights, etc.) as well as controlled active LED illumination near the fingers. Ambient correction using TCI (sec¬ 
ond column) removes the contribution of ambient illumination to the scene, yielding transmittance due solely 
to active LED illumination of which the spectral and power characteristics are known. See Media 1 for video 
results of ambient correction with varying illumination. 
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Figure 6: Example of the detrending procedure. A slow trend is estimated in the original signal, and subsequently 
removed to obtain a signal with a stable mean. This process removes variability in illumination due to movement 
and subtle changes in incident light. 
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Table 1: Experiment 1 statistical results for short-distance measurements. Pearson’s linear correlation coefficient 
(p) was computed between the PPGI and PPG signals for the whole 10 s signal for global accuracy, and for each 
heart beat instance (i.e., diastole-to-diastole) for local accuracy. A p -value was computed using a paired-sample 
t-test to test the hypothesis that there is a zero-mean difference between unprocessed and processed correlation. 
At the whole-signal level, the proposed system processing yielded greater than doubly improved accuracy over 
the unprocessed signal. At the local level, statistically significant results (a = 0.05) were achieved for each 
participant, indicating significant improvement with processing. 



pi 

P2 

P3 

P4 

P5 

Unprocessed, whole 

0.35 

0.31 

0.26 

0.14 

0.12 

Processed, whole 

0.94 

0.79 

0.82 

0.49 

0.56 

Unprocessed, blocked* 

0.72 ± 0.16 

0.51 ± 0.17 

0.65 ± 0.16 

0.22 ± 0.32 

0.28 ± 0.36 

Processed, blocked* 

0.96 ± 0.03 

0.79 ± 0.12 

0.85 ± 0.05 

0.51 ± 0.21 

0.59 ± 0.32 

p -value 

< 0.0001 (n = 13) 

< 0.0001 (n = 14) 

< 0.0001 (n = 9) 

0.040 (n = 10) 

0.020 (n = 13) 


* p ± a 


Table 2: Experiment 2 statistical results for long-distance measurements. Pearson’s linear correlation coefficient 
(p) was computed between the PPGI and PPG signals for the whole 10 s signal for global accuracy, and for each 
heart beat instance (i.e., diastole-to-diastole) for local accuracy. A p -value was computed using a paired-sample 
t-test to test the hypothesis that there is a zero-mean difference between unprocessed and processed correlation. 
At the whole-signal level, the proposed system processing yielded greater than doubly improved accuracy over 
the unprocessed signal. At the local level, statistically significant results (a = 0.05) were achieved for each 
participant except P3 (p = 0.056), indicating substantial improvement with processing. 



pi 

P2 

P3 

P4 

P5 

Unprocessed, whole 

0.18 

0.21 

0.14 

0.09 

0.10 

Process, whole 

0.62 

0.61 

0.45 

0.66 

0.39 

Unprocessed, blocked* 

0.38 ± 0.24 

0.44 ± 0.21 

0.13 ± 0.28 

0.32 ± 0.20 

0.12 ± 0.29 

Processed, blocked* 

0.72 ± 0.19 

0.63 ± 0.11 

0.43 ± 0.22 

0.70 ± 0.14 

0.33 ± 0.36 

p -value 

< 0.0001 (n = 15) 

0.0088 (n = 13) 

0.056 (n = 9) 

0.00040 (n = 10) 

0.0079 (n = 11) 


*/.i =t a 
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